AD-A125  902  LOCALLY  ONE  DIMENSIONAL  NUMERICAL  METHODS  FOR  1/1 

MULTI-DIMENSIONAL  FREE  SURFACE  PROBLEMS(U)  GEORGIA  INST 
OF  TECH  ATLANTA  G  H  MEYER  MAR  83  ARO- 16373 . 4-MA 
UNCLASSIFIED  DAAG29-79C-0145  F/G  12/1  NL 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  HiiHfAO  ft  >MM‘AKI>x  ;  i<. .  i 


OHC  FILE  COPY  JA  125902 


Unclass i f ied _ 

SECURITY  CLASSIFICATION  Qf  THIS  PAGE  (Wien  0«f»,8iil»«O _ _  .... _ 

_ REPORT  DOCUMENTATION  PAGE  BepoIePc^e?SgNfSorm 

I.  REPORT  NUMBER  "  |2.  GOVT  ACCESSION  MO.  J.  RECIPIENT1*  CAT ALOG  NUMBER 

16373. 4-ha _ A  h *4 /  'J- 5*  AllJz-  _ 

«.  TITLE  (m*d  Subtitle)  »•  TYEI  Of  MW*T  t  «*IOO  COVIMD 

Locally  One  Dimensional  Numerical  Methods  for  ,  .  iw  Ro 

Multi-Dimensional  Free  Surface  Problems 


t>. 4/  ^ £ 


Locally  One  Dimensional  Numerical  Methods  for 
Multi-Dimensional  Free  Surface  Problems 


t.  author**; 

Gunter  H.  Meyer 


•  .  CONTRACT  ON  GRANT  NUMBERS 

DAAG23  79  C  0145 


».  performing  organization  name  and  adores* 

Georgia  Institute  of  Technology 
Atlanta,  GA  30332 

«i.  controlling  office  name  and  address 
U.  S.  Army  Research  Office 
Post  Office  Box  12211 
Research  Triangle  Park,  NC  27709 

Is.  MONITORING  AGENCY  name  S  AOORCSSfff  dlllerent  I 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  *  WORK  UNIT  MUMPERS 


12.  REPORT  DATE 


IS.  NUMBER  OF  PACES 


Controlling  OtU ee)  |  IS.  SECURITY  CLASS,  (ol  1*1  •  report) 


I  is.  DISTRIBUTION  STATEMENT  (at  Olio  Keport) 


.i.  K  CLASSIFICATION.  OOVNGRADINO 
SCHEDULE 


Approved  for  public  release;  distribution  unlimited. 


1 17.  DISTRIBUTION  STATEMENT  (ol  tho  Mt IrKI  HlRM  In  Block  20,  M  dll  term!  Bern  Report) 


tvi;AA  .  S63  g 


ia.  supplementary  notes 

The  view,  opinions,  and/or  findings  contained  in  this  report  are  those  of  the 
author (s)  and  should  not  be  construed  as  an  official  Department  of  the  Army 
position,  policy,  or  decision,  unless  so  designated  by  other  documentation 

to.  KEY  WOROS  (Continue  an  warn  eldo  U  naeaaaarr  end  Identity  by  block  number) 

numerical  analysis 
free  boundary  problems 
partial  differential  equations 
algori thms 


Research  on  the  numerical  solution  of  free  boundary  problems  for  partial  different  al 
equations  with  locally  or  sequentially  one-dimensional  methods  has  been  supported 
by  the  U.S.  Army  Research  Office  through  two  consecutive  3"Y«ar  research  contracts 
The  time  and  resources  provided  have  made  it  possible  to  develop  a  reasonably 
comprehensive  mathematical  theory  and  flexible  numerical  algorithms  on  which  to  ba  e 
current  computational  methods  and  future  research.  During  the  first  3-year  period, 
the  method  of  fractional  steps  and  the  method  of  lines  were  applied  to  elliptic. 


00  1473  EDITION  OF  *  MOV  BS  IS 


070 


UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  RAM 


. 

CCv 

Dole  Metered) 


StCUWlTY  CLAttlFICATiQN  or  I  Hit  MOggw  PW  l»l— Q 


and  parabolic  free  boundary  problems.  During  the  past  three  years  work  was 
directed  toward  demonstrating  the  flexibility  of  the  method  of  lines  for 
increasingly  complex  problems,  on  examining  the  behavior  of  certain  ill-posed 
elliptic  free  boundary  problems,  and  on  establishing  a  mathematical  theory 


for  sequentially  one-dimensional  algorithms 


final  Report 

(TWENTY  COPIES  REQUIRED) 

1.  ARC  PROPOSAL  NUMBER: _ 

2.  PERIOD  COVERED  BY  REPORT:  8/1/79  -  12/31/82 _ 

3.  TITLE  OF  PROPOSAL:  Locally  one  dimensional  numerical  methods  for  multi-dimensional 

free  surface  problems 

4.  CONTRACT  OR  GRANT  NUMBER:  DAAG  29-79-C-0145 _ 

5.  NAME  OF  INSTITUTION:  Georgia  Institute  of  Technology _ 

6.  AUTHOR (S)  OF  REPORT:  Gunter  H.  Meyer _ 

7.  LIST  OF  MANUSCRIPTS  SUBMITTED  OR  PUBLISHED  UNDER  ARO  SPONSORSHIP 
DURING  THIS  PERIOD,  INCLUDING  JOURNAL  REFERENCES: 

Attached 


8.  SCIENTIFIC  PERSONNEL  SUPPORTED  BY  THIS  PROJECT  AND  DEGREES  AWARDED 
DURING  THIS  REPORTING  PERIOD: 


Final  Report  on  research  performed  under  Contract  DAAG29-79-C-0145  on 
"Locally  one  dimensional  numerical  methods  for  multi-dimensional 
free  surface  problem." 

Research  on  the  numerical  solution  of  free  boundary  problems  for  partial 
differential  equations  with  locally  or  sequentially  one-dimensional  methods 
has  been  supported  by  the  U.S.  Army  Research  Office  through  two  consecutive 
3-year  research  contracts.  The  time  and  resources  provided  have  made  it 
possible  to  develop  a  reasonably  comprehensive  mathematical  theory  and 
flexible  numerical  algorithms  on  which  to  base  current  computational  methods 
and  future  research. 

During  the  first  3-year  period,  the  method  of  fractional  steps  and 
the  method  of  lines  were  applied  to  elliptic  and  parabolic  free  boundary 
problems.  Both  methods  lead  to  sequentially  one-dimensional  algorithms  for 
multi-dimensional  partial  differential  equations.  They  are  judged  to  be 
conceptually  simple,  easily  coded,  and  little  affected  by  complex  free 
boundary  conditions.  A  more  detailed  account  of  this  phase  of  the  research 
is  contained  in  the  final  report  on  Grant  DAAG-29-76-G-0261  (copy  attached). 

During  the  past  three  years  work  was  directed  toward  demonstrating  the 
flexibility  of  the  method  of  lines  for  increasingly  complex  problems,  on 
examining  the  behavior  of  certain  ill-posed  elliptic  free  boundary  problems, 
and  on  establishing  a  mathematical  theory  for  sequentially  one-dimensional 
algorithms.  In  particular,  the  locally  one-dimensional  solver  was  applied 
to  coupled  heat  and  man  transport  during  alloy  solidification  [  4  ].  This 
work  is  preparatory  to  applying  locally  one-dimensional  methods  to  multi¬ 
dimensional  reaction-diffusion  systems.  With  the  same  goal  in  mind  monotone 
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iterations  were  applied  to  elliptic  free  boundary  problems  with  nonlinear 
source  terms.  In  this  setting  a  complete  convergence  analysis  can  be  given 
for  the  method  of  lines  when  applied  to  free  boundary  problems  which  are 
equivalent  to  variational  inequalities. 

At  the  same  time  Hele-Shaw  flow  in  a  porous  medium  was  examined.  The 
Hele-Shaw  injection  problem  is  readily  solvable  and  completely  analyzable 
with  variational  inequalities.  In  contrast,  the  suction  problem  is  ill-posed, 
has  no  variational  structure  and  has  not  yielded  to  analysis.  Our  numerical 
method,  on  the  other  hand,  remains  applicable  but  requires  careful  attention 
to  implementation  C  5  ].  An  attack  on  the  suction  problem  with  conformal 
mapping  [  6  ]  yielded  some  special  analytic  solutions  against  which  numerical 
results  can  be  compared. 

Among  our  current  projects  we  are  examining  whether  surface  tension 
effects  can  stabilize  the  suction  problem.  This  work  requires  extensive 
numerical  experiments  on  a  Cyber-type  computer.  At  the  same  time,  we 
have  begun  work  on  implementing  locally  one-dimensional  methods  on  a 
micro-computer. 
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Gleichungen,  Seminar  uber  die  Linienmethode  und  Anwendunzen, 

K.  J.  Schwenzfeger,  edt.,  Hochschule  der  Bundeswehr  Munich,  1982. 
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GlINTER  H.  MEYER 


Final  report  on  work  performed  under  Grant  DAAG  29-76-G-0261  on 
"Locally  one  dimensional  numerical  methods  for  multidimensional 

free  surface  problems” 

The  applicability  of  locally  one  dimensional  methods  to  multidimen¬ 
sional  free  boundary  problems  was  examined.  Such  problems  require 
the  solution  of  an  elliptic  or  parabolic  field  equation  over  a 
domain  which  is  unknown  a  priori  and  which  must  be  determined  from 
over-prescribed  boundary  data.  Our  main  objective  was  to  identify 
a  numerical  method  for  such  problems  which  is  widely  and  easily 
applicable,  but  which  at  the  same  time  has  a  solid  mathematical 
foundation. 

Initially,  we  examined  the  method  of  fractional  steps  which  for  a 
multidimensional  diffusion  equation  leads  to  a  sequence  of  loosely 
coupled  one  dimensional  free  boundary  problems  in  alternating  ortho¬ 
gonal  directions.  In  a  report  prepared  just  prior  to  the  grant 
period  this  method  is  introduced  and  used  together  with  invariant 
imbedding  to  solve  a  variety  of  Stefan  type  free  boundary  problems. 
The  method  works  well  and  allows  comparable  numerical  resolution 
in  all  space  directions.  However,  it  does  impose  certain  restric¬ 
tions  on  the  shape  of  the  free  boundary  which  serverely  limit  the 
applicability  of  the  ADI  method. 

To  overcome  this  limitation  we  began  working  with  the  method  of 
lines  based  on  a  simple  finite  difference  discretization  of  time  and 
of  all  but  one  space  variable.  The  original  multidimensional 
equation  is  thus  replaced  by  a  system  of  one  dimensional  problems 
along  parallel  lines.  This  system  is  solved  in  line  SOR  fashion 
so  that  again  a  sequence  of  uncoupled  one  dimensional  problems  re¬ 
sults.  As  above  we  apply  invariant  imbedding  to  each  one  dimen¬ 
sional  problem.  A  description  of  this  approach  appears  in  II).  A 
comparison  with  the  fractional  step  method  for  a  two  dimensional 
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ablation  problem  [2]  shows  that  only  relatively  few  lines  compared 
to  the  number  of  such  points  along  each  line  should  be  used  and 
hence  that  the  numerical  resolution  is  not  equally  fine  in  all 
directions.  Nonetheless,  even  in  this  severe  test  problem  the 
method  of  lines  performs  as  well  as  the  ADI  method.  Because  free 
surfaces  of  more  general  shape  are  admissible  further  work  has  con¬ 
centrated  on  the  method  of  lines/SOR/invariant  imbedding. 

The  success  of  this  method,  even  in  the  crude  implementation  em¬ 
ployed  as  far,  was  good  motivation  to  establish  the  mathematical 
validity  of  the  algorithm.  As  a  first  step  convergence  of  the  it¬ 
erative  numerical  method  and  the  discretization  was  established 
for  Poisson's  equation  subject  to  nonlinear  boundary  data  on  a 
fixed  domain  [3].  Recently,  the  arguments  were  extended  to  estab¬ 
lish  convergence  for  an  elliptic  free  boundary  problem  arising  in 
hydrodynamic  lubrication  [4] .  As  far  as  we  are  aware  this  is  the 
first  time  that  a  multidimensional  front  tracking  algorithm  can  be 
completely  analyzed.  Current  work  is  concerned  with  convergence 
for  the  equations  arising  in  fluid  flow  through  a  porous  dam.  This 
problem  is  characterized  by  considerably  more  complicated  boundary 
conditions  than  present  in  the  lubrication  problem. 

Although  we  are  primarily  concerned  with  multidimensional  problem 
the  one  dimensional  solver  is  at  the  heart  of  our  method  and  we 
continue  to  gather  experience  with  the  invariant  imbedding  method. 
An  application  to  a  problem  with  distributed  boundary  conditions 
is  described  in  [7} . 

During  the  grant  period  we  particapted  in  two  conferences  dedicated 
to  free  and  moving  boundary  problems.  A  survey  over  numerical  meth¬ 
ods  for  Stefan  problems  and  a  comparison  with  the  method  of  lines 
was  presented  in  [5],  while  an  application  of  our  methods  to  a 
variety  of  practical  problems  including  heat  and  mass  transfer, 
Hele-Shaw  flow,  and  three  dimensional  ablation  is  described  in  (6], 
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Two  additional  publications  were  prepared  for  the  1977  and  1978 
Army  Numerical  Analysis  Conferences.  A  revision  of  the  first  re¬ 
port  was  published  as  [2]  while  the  material  of  the  second  is  con¬ 
tained  in  [6]. 

Colloquium  talks  on  our  work  were  presented  during  the  past  three 
years  of  several  universities  and  laboratories  in  the  U.S.,  Canada, 
Austria,  Germany  and  Great  Britain. 

In  summary,  the  basic  premise  that  sequential  one  dimensional 
methods  provide  a  useful  tool  for  multidimensional  free  boundary 
problems  has  proved  to  be  correct.  It  remains  to  establish  firmly 
the  limits  of  their  applicability  and  to  provide  an  analysis  for 
a  wider  class  of  problems  than  was  possible  so  far.  First  and 
foremost,  however,  the  software  aspects  of  our  methods  deserve 
attention  in  order  to  improve  the  efficiency  of  our  methods. 


During  the  grant  period  the  following  graduate  student  assisted  in 
our  project. 
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F.  Gerwig,  M.S.  from  Georgia  Tech,  1977 
D.  Dvorak,  graduate  student  in  mathematics 
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